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We present a statistical description of the propagation of short pulses in long 
optical fibers, taking into account the Kerr and nonlocal nonlinearities on 
an equal footing. We use the Wigner approach on the modified nonlinear 
Schrodinger equation to obtain a wave kinetic equation and a nonlinear dis- 
persion relation. The latter exhibit that the optical pulse decoherence reduces 
the growth rate of the modulational instability, and thereby contribute to the 
nonlinear stability of the pulses in long optical fibers. It is also found that the 
interaction between spectral broadening and nonlocality tends to extend the 
instability region. 
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The nonlinear propagation of pulses in optical fibers has attracted a great deal of interest 
since the early seventies^, and research into the application and theory of this field is still 
growing-. The evolution of coherent weakly nonlinear optical pulse envelopes is given by 
the cubic nonlinear Schrodinger equation (CNLSE)^ involving the Kerr nonlinearity. The 
CNLSE admits bright, dark and gray solitons, which are used for ultrahigh-speed optical 
communications without pulse broadening and spectral dilution 4 -. Even so, it has become 
clear that the effects of noise in fibers and amplifiers could alter the pulse properties in 
significant ways. Thus, it is of crucial importance to obtain qualitative and quantitative 
estimates of the effects of different types of incoherent perturbations 5,6 on the optical pulse 
propagation. Recently, the Wigner transform technique^ in nonlinear dispersive media has 
been used to analyze Landau-like damping*^, partially coherent higher order dispersive 
effects^, the importance of the incoherence spectrum^, and the influence of incoherence on 
the modulational instabilit y 13 ^ 4 for cases involving the cubic Kerr nonlinearity. 

However, there are other important nonlinearitie a 15 i 16 i 17 (e.g. saturation and higher order 
nonlocal nonlinearities) which can compete with the cubic nonlinearity in optical fibers. The 
combined influence of the cubic and nonlocal nonlinearities on the modulational instability of 
a constant amplitude coherent optical pulse has been examined by Shukla and Rasmussenii. 
In this Letter, we present a statistical description of partially incoherent pulses in long optical 
fibers, taking into account the Kerr and nonlocal nonlinearities on an equal footing. We use 
the Wigner approach and deduce a wave kinetic equation from which a nonlinear dispersion 
relation (NDR) has been derived. The NDR is then analyzed to demonstrate the effect of 
random noise on the modulational instability of incoherent optical pulses. It is found that 
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the optical pulse decoherency can contribute to the nonlinear stability of pulses in optical 
fibers. 

Given the electric field E(z, t) exp(ik z — iu t) of the optical pulses, the evolution the 
pulse envelope E in the slowly varying envelope limit, i.e. k ^> (d z — 2k' d t ), is governed 
by 17 

i(d z + T)E + ad 2 t E + (3IE + i^d t (IE) = 0, (1) 

where we have introduced the parameters a = —k'l/2, (3 = 712/co/no, 7 = 2n 2 /c, and 
r = koXo/ n o- Moreover, the prime denotes differentiation with respect to u>o, the intensity 
parameter is given by I = \E\ 2 , the refractive index is n(uj ,I) = n + ixo + n 2 I, n Q = n(uo), 
and and Xo = x( w o) represents losses in the medium. 

In order to take the effects of partial coherence into account, we define the space-time 
correlation function for the electric field as C(z+, Z-, t + , t_) = E*(z+,t + )E(z-,t-), where 
z± = z ± C/2 and t± = t ± r/2. Then, the Wigner distribution function of the optical pulse 
is given by*^ 

F{z,t,k,u) = j±p J d(dTe^-^C(z + ,z„,t + ,t„), (2) 



such that 



1 
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I(z,t) = ( 2 I dkdu F(z,t,k,u). (3) 



Thus, from Eq. (Q) the evolution equation for the Wigner function @ corresponding to 



the envelope field E becomes (see also Ref. 



18) 



2uad t F - d z F + 2(31 sin (j F + 7 \-d t [/cos (§ S£) F~\ + 2^/sin (§ d t dj) f} = 2rF(4) 

where we have performed the Wigner transformation over the time domain. Here the arrows 



denotes direction of operation, and the operator functions are defined in terms of their 
respective Taylor expansion. The system of equations Q and (J1J) determines the evolution 
of short partially coherent optical pulses in nonlinear media. 

In order to analyse the modulational instability and the effects of the terms due to a 
nonzero 7, we make the ansatz F(z, t, u) = F (uj) + Fi(uj) exp(iKz — iQt) + c.c, where c.c. 
denotes the complex conjugate, and \Fx\ <C F . Moreover, since we are interested in the 
short-pulse effects, we will for simplicity neglect the loss term T in Eq. (@J, in order to obtain 
clearly interpretable results. 1 Expanding Eq. (j3J) in terms of this ansatz, and using Eq. @, 
we obtain 

1 f , [P + jju + n/2)] Fp(u - n/2) -[(3 + j(fj - n/2)] F (uj + Q/2) 
2anJ lu + (K- 7 QI )/2aQ ' 1 j 

where Io = f duF (uj). Equation (jSJ) represents the NDR for a short optical pulse, where 
the pulse may have spectral broadening and partial coherence. 

In the case of a mono-energetic pulse, we have F (u) = I 5(ui — Qo), where O corresponds 
to a frequency shift of the background plane wave solution, and the NDR © gives 17 

K = 2(7/0 - «fi )^ ± [7 2/ o fi2 + « 2fi4 - 2a W + 7^o)^ 2 ] 1/2 • (6) 

In practice however, the wave envelope will always suffer perturbations due to various 
noise sources, e.g. fiber and amplifier noise. A noisy environment may cause the pulse field 
to attain a random component in its phase. Thus, if the phase <p(x) of the electric field 
varies stochastically, such that the ensemble average of the phase satisfie a 20 i 21 (exp[— i(f(t + 

1 It should be stressed that in certain applications, the losses may not be small, and the T term should 
under these circumstances be kept. As noted by Shukla and Rasmussenii the effect of the loss term is to 
damp the pulse according to exp(— 2Tz) as it propagates through the fiber. 
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t/2)] exp[i(p(t — t/2)}) = exp(— f2r|r|), the background Wigner distribution is given by the 
Lorentzian spectrum 

IT (iU — i] ) + "T 

where corresponds to the width of the spectrum. Then, the NDR (j3J) takes the form 



27 [if — 7 / ^ + an(n - in T )} - 2 a pn 

which has the solution 

if = 2 [7/0 - a(n - itor)] n ± [7 2 / 2 ^ 2 + a 2 n 4 - 2al ((3 + 7(^0 - itt T ))n 2 ] 1/2 • (9) 

This solution generalizes the result © to the case of a random phase background envelope 
field. Equation © clearly shows that the width gives a nontrivial contribution to the NDR. 
We note that when 7 = 0, we may define the growth rate k according to K = — 2aQ Q — in, 
and the width Qt then gives rise to a Landau like damping from Eq. (jHJ). 

When 7 is non-zero, the growth/damping behavior becomes considerably more complex, 
with new instability regions. Letting / = 7 2 Jq + a 2 Q 2 — 2al (j3 + 7^0), and assuming 
<C f/ a jla, we obtain the approximate expression 

K/VL « 27/ - 2att ± f 1/2 + 2iatt T ± ia^I VL T / f 1/2 (10) 

from Eq. ©. When a > 0, and 2al ((3 + 7^0) > 7 2 /q + « 2 ^ 2 , we have / < 0. Denoting 
the growth rate by k = — Im(if), we obtain k = \f\ 1 ^ 2 — 2iaVtx- Thus, as expected, the 
coherence spread gives rise to a smaller growth rate for the modulational instability of 
incoherent optical pulses. We note that this instability occurs also when 7 = 0. 

On the other hand, if a > 0, but 2al ((3 + 7^0) < 1 2 Iq + « 2 ^ 2 , or a < 0, so that / > 0, 
a new novel effect is present due to a nonzero 7. We have k = (7/0/ f 1 ^ 2 — 2) o&It- Thus, 

5 



a short pulse in conjunction with a finite statistical spread f2y could give rise to a shift in 
the damping due to the decoherency of the pulse, which hence implies a shift also in the 
growth rate. This effect can be seen in Fig. [TJ where we have plotted n as given by the full 
dispersion relation Q for the frequency shift fl = 0. We have used the rescaling I — > (3I , 
$It — > y/a£l T , n — > y/aCl, and 7 — > j/(P^/a) = V2 n /(ck \k'Q I 1 / 2 ). We note that not only 
is the damping shifted, but the instability regions is also extended, and quite significantly for 
higher values of j/(/3y/a). Since r y/(/3y/a) oc \D\~ 1 ^ 2 , where D is the dispersion parameter 
commonly used in fiber-optics, the value of the normalized non-locality strength may become 
large, as D can be designed to be very close to zero for certain wavelengths 22 . Thus, the 
novel coupling between spectral broadening and nonlocality should be possible to measure 
using a suitable setup. 

To summarize, we have presented an investigation of the modulational instability of inco- 
herent optical pulses in a nonlinear optical medium that contains the Kerr and higher order 
nonlocal nonlinearities on an equal footing. By using the Wigner transform, we have derived 
a wave kinetic equation for incoherent pulses from the generalized nonlinear Schrodinger 
equation. The wave kinetic equation is further exploited to obtain a nonlinear dispersion 
relation, which exhibits new features of the modulational instability. We find that the deco- 
herence of the optical pulses reduce the modulational instability growth rate due to a spatial 
damping caused by the broad optical pulse spectrum. However, the combined effect of a 
random phase and a non-local nonlinearity is to extend the instability region as compared 
to the case of a monochromatic spectrum. Thus, the present result thus contribute to the 
nonlinear stability of incoherent optical pulses in long optical fibers. 
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Fig. 1. The effects of spectral broadening and non-locality. Normalizing the variables K 
and Q, as well as the parameters Qt and 7, such that a — (3 — 1, we have plotted the 
imaginary part k = — lm(K) as a function of Q, when the frequency shift f2 is put to 
zero. From the peaks of the curves downwards, we have used I = 0.5, and the full curve 
represents Qt = 7 = 0, and shows the regular modulational instability growth rate. The 
next curve (dashed) gives k for fl T = 0.1 and 7 = 0, while the third (dashed-dotted) curve 
uses f2 T = and 7 = 1, and the fourth (dashed-dotted) curve has f2 T = 0.1 and 7 = 1. The 
last two curves (dashed and dotted, respectively), where Qt = 0.1, 7 = 1.9, and f2 T = 0, 
7 = 1.9, respectively, clearly shows the character of the combined effect of broadening and 
non-locality, namely a widening instability region. 
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